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Abstract. Excitonic transitions offer a possible route to ultrafast optical spin 
manipulation in coupled nanostructures. We perform here a detailed study of the three 
principal exciton-mediated decoherence channels for optically-controlled electron spin 
qubits in coupled quantum dots: radiative decay of the excitonic state, exciton-phonon 
interactions, and Landau-Zener transitions between laser-dressed states. We consider a 
scheme to produce an entangling controUed-phase gate on a pair of coupled spins which, 
in its simplest dynamic form, renders the system subject to fast decoherence rates 
associated with exciton creation during the gating operation. In contrast, we show that 
an adiabatic approach employing off-resonant laser excitation allows us to suppress all 
sources of decoherence simultaneously, significantly increasing the fidelity of operations 
at only a relatively small gating time cost. We find that controUed-phase gates accurate 
to one part in 10^ can realistically be achieved with the adiabatic approach, whereas 
the conventional dynamic approach does not appear to support a fidelity suitable for 
scalable quantum computation. Our predictions could be demonstrated experimentally 
in the near future. 
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1. Introduction 

Quantum computation (QC) promises a significant speedup for certain classes of 
problems [H [2l [3] as well as an efficient way of simulating quantum systems [1]. The 
broad range of expertise and knowledge that exists in fabricating and characterising 
semiconductor structures has naturally led to a number of proposals for quantum dot 
(QD) QC being put forward. The spin of a single electron confined within a QD 
has been suggested as a scalable qubit [5], while another possibility is to define the 
computational basis states as the presence or the absence of a confined electron-hole 
pair (an exciton). While excitonic qubits can be conveniently optically addressed, they 
also interact strongly with their surroundings and thus suffer from rapid decoherence 
(typically within a nanosecond [6]). Spin-based qubits, on the other hand, have much 
longer coherence times (up to a millisecond [7j) but can suffer from the fact that inter- 
spin interactions tend to be weak. Hybrid schemes have therefore emerged that propose 
to store the qubit in a long-lived localized electronic spin state, moving selectively to an 
excited state spin-exciton (trion) representation only when a gate is actually performed. 
Such schemes hope to marry the advantages of the two different representations [HIIH], 
though it should be borne in mind that any excitonic component of the system state 
will still suffer from rapid decoherence, and hence significant excited state population is 
ideally avoided. 

For universal QC it suffices to have arbitrary single qubit operations and at least 
one entangling two qubit gate [1]. Single qubit operations are simply rotations of 
the state vector on the Bloch sphere and proposals based both on stimulated Raman 
adiabatic passage and direct optical excitation exist for spins in QDs [lOl [11], [T2] . 
Decoherence processes have been thoroughly studied in these cases [HI [131 El [15] . 
leading to the prediction that high-fidelity single-qubit operations should be feasible. In 
many physical systems, the easiest entangling operation to implement is a controUed- 
phase (CPHASE) gate that leaves three of the four possible two-qubit computational 
basis states unaltered, while generating a phase of vr on the fourth (e.g. |00) — |00), 
|01) |01), |10) — s> |10), |11) — 111)) [I]- Various schemes to realise such a gate exist 
for spins in coupled QDs [H [161 [HI; of which adiabatic optical control is a particularly 
promising example as it naturally suppresses excited state population throughout the 
operation [8], [18]. The technique relies on varying external control parameters of the 
system, here the intensity and frequency of a laser, sufficiently slowly compared to 
the characteristic timescales of the system itself that the entire population remains in 
instantaneous eigenstates throughout the gating operation. This approach overcomes 
light-heavy-hole-mixing in the QD valence band [HI [E] and is also thought to be robust 
to phonon-induced dephasing mechanisms [8]. While the adiabatic approach seems 
extremely promising as perhaps the best way to achieve rapid control of solid state 
spins, it does have potential issues like every approach. It is therefore essential to 
perform a comprehensive study of all major noise sources during a two-qubit gate. 
Here we provide such a study, and conclude that the approach is in fact remarkably 
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robust. We address the principal decoherence channels of radiative recombination and 
carrier-phonon interactions by deriving a Markovian master equation (ME) for the laser- 
driven two-dot system, accounting for both processes. We will use the ME to assess 
the performance of the adiabatic approach within experimentally accessible parameter 
regimes. We also derive conditions to avoid non-adiabatic Landau-Zener (LZ) transitions 
between eigenstates [191 1201 [211 [22], an effect that poses an additional error source. In 
order to put the results in context, we will consider a simple dynamic approach based 
on resonant Rabi flopping as a comparison. 

2. The System 

By observing Rabi oscillations, optical coherent control of excitons in single QDs has 
been demonstrated by a number of groups in recent years [23l [24l [25]. Further, 
excitonic interactions between QDs have been observed experimentally [261 [271 [28] and 
optical conditional logic between two excitonic qubits within a single QD has been 
accomplished [29]. QD spin qubit Rabi flopping has been achieved both through direct 
spin resonance (on the sub-microsecond time scale) [30] and optically (on the picosecond 
time scale) [31]. Single-shot read-out by spin to charge conversion [32] has also been 
achieved. In addition, Greilich et al. [33] have performed optical measurements of spin 
coherence and Xu et al [31] have done fast spin cooling using excitons. 

We consider a realistic state-of-the-art system: Two adjacent self- assembled QDs (I 
and II respectively) with a distinct heavy- and light-hole valence band structure, such as 
those grown on a GaAs substrate [35]. The dots could be placed side-by-side, though are 
more usually grown in stacks (see Ref. [36] for example). We consider QDs that are small 
and have a strong confinement potential that dominates over any intra-dot Coulomb 
interactions. They are each doped such that a single excess electron, the spin of which 
embodies the qubit, permanently occupies the lowest energy state of the conduction 
band. The qubit basis is defined as |0) = | |) = | — 1/2) and |1) = | |) = |l/2), 
where jm^) is the spin projection, with = ±1/2. The z direction also defines the 
quantization axis for the hole spins (see below). Our qubit representation is different 
to that used by Imamoglu et al. [9], where an applied magnetic field defines the qubit 
basis in the x-direction. Further, since our scheme relies on direct inter-dot interactions 
rather than the cavity mediated Raman transitions needed in Ref. [9], we require only 
a single off-resonant laser pulse to effect a CPHASE gate. 

In order to couple the two spin qubits we exploit the fact that trion-trion inter-dot 
Coulomb interactions are expected to be relatively strong [371 [IS] such that we can 
address the system optically through spin-selective exciton creation [H [371 [iSj- This 
effect, which has recently been observed experimentally [3H1 [39], relies on the Pauli 
exclusion principle for the lowest energy conduction band electron state (see Fig. [1]). 
This state can be occupied by at most two electrons, each of opposite spin orientation, 
with the electron qubit permanently occupying one of the slots. In the valence band, 
the heavy and light hole energy levels split due to their differing effective masses. Heavy 
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Figure 1. Pauli blocking effect: cr+ light cannot excite electron hole pairs if the qubit 
is in the state |0) with = — 1/2 (left hand side). On the other hand, if the qubit 
is in the |1) configuration the excitation is possible, leading to a three particle trion 
state \X) (right hand side). 



holes occupy the lowest energy states with spin = ±3/2 (whereas the light holes have 
spin = ±1/2) and are the only valence band states considered here (see Ref. [18] 
for detailed calculations of the effects of light-heavy hole mixing). A a"*" circularly 
polarized laser incident on the QD carries photons with angular momentum / = ±1. 
Thus, if the qubit is in state |0), with spin = —1/2, promoting another electron to 
the conduction band is incompatible with the conservation of angular momentum and 
no transition can occur. Conversely, for a qubit in state |1) (m^ = 1/2) the transition is 
allowed and leads to the creation of a three-particle trion state \X) = \ ], where 
the arrows symbolize the spin projection of the electron, and the triangle of the hole 
(|V) = I — 3/2)), respectively. 

That Coulomb interactions between trions on the two adjacent dots mediate a 
coupling of the two spin qubits can be seen by considering the full system Hamiltonian. 
We assume that the qubit states are degenerate and write the Hamiltonian of the two 
QDs and a single classical laser field in the basis {|0), |1), \X)} {h = 1) as 

H{t) = uJo{\X){X\ ® 1 ± 1 ® \X){X\) 

+ncosuJit{\l){X\ 01 + 10 \1){X\ + h.c.) 

+Vxx\XX){XX\ + Vf{\1X){X1\ + h.c). (1) 

where h.c. denotes the Hermitian conjugate, 1 is the identity operator, and uq is the 
exciton creation energy (assumed the same for both dots). Here, Q represents the 
coupling between the QD transition and laser mode and ui is the laser frequency, both 
of which may be time-dependent quantities. The dots are Coulomb coupled both by 
virtual photon exchange (Forster interaction) of strength Vp [ID], and by a biexcitonic 
dipolar coupling Vxx [37] . We consider typical inter-dot separations of 5 - 10 nm, close 
enough for both Vp and Vxx to be on milli-electronvolt scale. We also assume that we 
are in a regime where tunneling processes are suppressed |41j . 

To lowest nonzero order, the Forster coupling element is equivalent to the 
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interaction of two point dipoles situated on dots I and II [12], 

Vf = ^, ((ri) ■ (rn) - |j ((ri) ■ R) ((m) • R)) , (2) 

where R is the vector connecting the centre of the two dots, (ri/ii) is the position 
operator between the electron and hole on either dot I/II respectively, is the dielectric 
constant, and C the usual Coulomb term C = e^/47reo. The sign of Vp is determined by 
the relative orientation of (rj) and (rn) with respect to R. Ultimately, the orientation 
of (rj/ji) depends on the lattice vector structure jl2] and is hard to predict. However, 
for two identical dots in the same crystal structure, (ri) and (rn) should be parallel 
and have the same magnitude, such that their joint orientation to R is the sole factor 
determining whether Vp is positive or negative. While the sign of Vp is an immutable 
property of any given two-dot system, both positive and negative values of Vp should 
be possible in principle by a suitable geometrical arrangement. 

It can be seen from Eq. ([1]) that the Hamiltonian decouples into four non- 
interacting subspaces [16]: Ho = {|00)}, Hi = {|01), |0X)}, Hr = {|10), |X0)} and 
^2 = {|11), \IX), \X1), \XX)}; this is a direct result of the Pauli-blocking effect. To 
implement the CPHASE operation, we need to achieve a net phase shift of tt on the 
input state |11), and are therefore primarily interested in the dynamics of 7^2- When 
we do need to consider dynamics in the two level subspaces Tii and T-Ci> we can appeal 
to the results of Ref. [15], which analyzed an isomorphic Hamiltonian structure. 

Proceeding as in [16], HS], we first transform 0.2 into the basis of its eigenstates in the 
absence of driving {n = 0): |11), |V'+) = {\IX) + \X1))/V2, |z/;_) = (|1X) - \X1))/V2 
and \XX), giving 

H2{t) = {uo - Vp)\ij^){ij^\ + {uo + Vp)\ij+){i,+ \ 

+V2ncosuit{\ll){ij+\ + \ij+){XX\ + h.c.) 

+ {2uo + Vxx)\XX){XX\. (3) 

We now move to a frame rotating with the laser frequency ui, detuned from the 
|11) ^ transition by an amount A = uq + Vp — uoi. Within the rotating wave 
approximation (RWA), Eq. Q becomes 

i/; = (A - 2Vp)\ij.){il)4 + A|^+)(^+| + (2A - 2Vp + Vxx)\XX){XX\ 
+ ^(|ll)(^<A+| + |^+)(XX| + h.c.). (4) 

Ref. [16] describes how to perform a dynamic CPHASE operation by applying a resonant 
2tt laser pulse to the |11) ^ transition (i.e. with A = 0), generating the required vr 
phase shift on state |11). In this case, certain conditions on the driving strength Vt must 
be satisified in order to suppress unwanted transitions to \XX) and transitions within 
Til and T-ty |T6]. However, the gate may also be operated adiabatically without those 
same constraints by slowly switching an off-resonant laser beam, forcing the system to 
follow its instantaneous eigenstates [IB]. The energy difference between states |11) and 
\'ipjf-) then allows for a phase accumulation on the |11) state, as can be seen by analysing 
the eigenstates in more detail. 
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Figure 2. Dressed states of subsystem 7^2- Parameters are: Vxx ~ 5 meV, 
Vp = 0.85 mcV, and = 1 meV. The inset sliows tire states and around tiie 
origin. Also shown are the dressed states obtained if \XX) is neglected. For positive 
A, these approximate and well. 



The three levels |11), and \XX) are coupled by the laser in Eq. (jlj). Moving to 
the diagonal basis results in three dressed states, which we shall label {|C-), |C-i-)) ICx)}, 
each of which are superpositions of the three bare basis states weighted according 
to certain mixing angles that are complicated functions of the system parameters. 
Nevertheless we can greatly simplify matters by an appropriate choice of parameters, 
and work with an approximate 2LS that describes the full dynamics with a high degree 
of accuracy. As shown in the inset of Fig. O for positive Vp = 0.85 meV, Vxx = 5 meV, 
and positive A, the energies of states and are very well approximated by taking 

|C_) ~ cose|ll) -sine|7/;+) 

IC+) ~ sine|ll) + cose|V'+), (5) 

where G = (1/2) arctan(-\/2^^/A). The agreement improves even further for larger A, 
which we shall demonstrate is desirable for the adiabatic scheme, or smaller driving fio- 
A negative Forster coupling Vp < simply shifts the dark state upwards in Fig. [2j 
At the same time, the avoided crossing between |(^+) and \(x) moves further to the right 
and hence interactions with the biexcitonic level are suppressed even more effectively. 
Therefore, we can safely assume that the perturbation caused by \XX) is negligible for 
a wide range of parameters. 

Let us now consider shining a chirped laser on the QD with temporal evolution 
A(t) and fl(t) that slowly and continuously changes 9 from to some value Qmax, and 
back to again. Since |11) and coincide for 6 = 0, any population in state |11) will 
adiabatically follow the instantaneous eigenstate of \(-) and return to |11) at the end of 
the pulse. As a consequence of its time evolution, the phase of an eigenstate changes 
according to \fi{t)) = exp[—iE^t)\fi{0)) , where is the relevant eigenenergy. Hence, 



Robust adiabatic approach to optical spin entangling in coupled quantum dots 



7 



at the end of the pulse, |11) picks up a phase of exp(— i^n) relative to |00), where 

T T 

(j)u = J dtE^_{t) = J dt^ [A(t) - v/A2(t) + 2f]2(t)] , (6) 



to a good approximation. Here, T is the pulse duration and Ec__{t) the energy of \QJ) 
relative to |00). Of course, states |01) and |10) also accumulate phase relative to |00) 
due to the laser-induced coupling to |0X) and |X0), giving 

T 



!^01,10 



1 

dt- 
2 L 



A'(t)- v/A'2(t) + fi2(i) ^ (7) 







respectively, where A'(t) = A(t) — Vp is the detuning in the Tii subspace. To achieve a 
gate locally equivalent to the CPHASE, ^oo — 0oi — 0io + 0ii = tt following the operation. 
This constrains the control parameters A(t), and T. 

For adiabatic following to hold, we must consider two further points. First, must 
change sufficiently slowly that LZ transitions between the dressed states do not occur; 
we will consider these in Section |5l Second, we must start and end with 6 = 0. In 
principal, both the laser detuning A(t) and the coupling strength VL{t) can be time- 
dependent, though it is essential to start with an off- resonant pulse and to slowly 
switch the coupling strength. It is then possible, although by no means necessary, 
to chirp the pulse and tune into resonance [8]. In this case the approach to resonance 
dominates the dynamics and the gating time of the operation. For the sake of simplicity 
and practicality, we will only consider here laser pulses with a constant detuning and 
Gaussian pulse envelopes. The duration of the pulse is then described by the parameter 
r: 

A(t) = A (8) 

=fioe-(*/")', (9) 

in contrast to previous work [HI [18] where chirped laser pulses were studied. Having 
described, in some detail, the internal dynamics of the driven system we now turn to 
the impact of the external environment on the CPHASE operation. 



3. Radiative decay 

The lifetime of excitons in a QD can be as long as a nanosecond [7], while coherent 
control should be possible on the picosecond timescale [Bj. However, in order to avoid 
phonon induced pure dephasing [15] and to maintain adiabaticity it can be advantageous 
to perform operations much more slowly. In this case the finite exciton lifetime becomes 
an important source of decoherence which cannot be neglected. 

Any population in an excited state of the system is susceptible to radiative decay 
at the rate of the inverse natural lifetime F, due to coupling with the vacuum radiation 
field. As discussed above, we need only consider the behaviour of an effective 2LS to 
capture such effects in 7^2, which we write i?2efr = ^H'+){'^+\ (in the basis {|11), \'ip+)}). 
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Figure 3. Radiative decay from manifold Ai{N) to /A{N — 1). Left: the micoupled 
basis. The energy difference between the two states in each manifold is the detuning 
A, the solid arrows correspond to absorption and stimulated emission processes, and 
the wavy arrows denote spontaneous emission. Right: Allowed spontaneous emission 
transitions between the dressed states. The energetic splitting in each manifold is the 
effective Rabi frequency ft' = y/ A'^ + 2^1'^ and the spacing between adjacent manifolds 
is the laser frequency oji. The leftmost emission process is at frequency uji + fl' , the 
two centre-lines emit at the frequency of the laser and the rightmost emission is at 
u!i — n': this gives rise to the famous MoUow triplet [47] . 



The system is coupled strongly to a laser mode and weakly to a reservoir of empty 
modes (the radiation field). In contrast to Section [2], we treat the laser mode quantum 
mechanically, assuming it to be monochromatic and in a coherent state with mean 
photon number (A^) ^ AA^ ^ 1. We also demand that the number of photons emitted 
into new modes as a result of fluorescence is much smaller than the width of the photon 
distribution AA^ ^ FT, with T being a characteristic dot-light interaction time. Under 
these assumptions, we can write the coupling term simply as r2(|ll)(^/'+| -|-h.c.) / a/2 |46j . 
This form agrees with the coupling term for a classical laser (within the RWA), and so 
we formally retrieve the relevant terms in Eq. (jlj). 

In absence of either the dot-laser coupling term or an environment, the energies of 
the combined QD and laser mode states {|11, A^ + 1), \ip+,N)} differ by the detuning 
A = ujq — uJi + Vp. If the detuning is small (A ^ uo), the states are energetically close 
to each other and can be conveniently grouped into a manifold Ai{N). Likewise, we 
group the pairs {|11, A^), \ij+,N - 1)}, {|11, A^ + 2), |^+, A^ + 1)} etc. into a ladder of 
manifolds differing in energy by steps ui. Introducing the dot-laser coupling mixes the 
bare eigenstates, again defining a dressed basis 

|_)^= cose|ll,A^+l) -sine|?/;+,Ar), (10) 

\+)n= sine|ll,A^+ 1) + cose|V^+,A^), (11) 

with 6 = (1/2) arctan (-\/2r2/A) as before, and we use the subscript A^ to denote the 
manifold Ai{N). Recall the similar set of eigenstates |— ) = cosO|ll) — sin Q\iIj^), 
1+) = sin ©I 11) + cosQ\ip^) that would be obtained by considering a classical laser. In 
contrast to the dressed states of Eqs. ffTOl) and ffTTl) . the classical treatment does not 
give any information on the number of photons in the laser mode nor on the manifold 
A^(A^). As a result, it appears at first glance that |— ) is truly the system's ground-state 
from which no relaxation is possible. Yet, it is clear that for 7^ both dressed states 
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have some excitonic character and a transition from each into either state of the adjacent 
manifold below, A4 [N— 1) , is possible. During this process a photon from the laser mode 
is emitted into a different, previously empty mode, by means of spontaneous emission 
(see Fig. [3]). As is easily seen in the dressed picture, three transitions are possible, 
emitting photons at frequencies ui, and lji ± Q', where the effective Rabi frequency 
= + 2f2^ is the energy separation between the dressed states in each manifold. 
In the context of the CPHASE gate, this implies that following the state ~ |— ) 
adiabatically does not leave the operation immune to radiative decay. On the contrary, 
if the pulse is tuned all the way to resonance, 6 = 7r/4, half the population will be 
susceptible to recombination, analogous in fact to the case of resonant excitation where 
half the population remains in state |11) and half in state on average. 

To estimate the impact of spontaneous emission on the gate fidelity we need to know 
the transition rates between the dressed states. Since our proposal relies on adiabatic 
following of the |— ) state, only its decay rate r„ needs to be considered. It is given by 
the sum of two processes: |— )Ar — > |— )Ar_i and \—)n |+)Af-i- The total decay rate is 
then 

r„ = V2To sin^ e cos^ 9 + v^Fq sin^ 6 = v^Fq sin^ 6, (12) 

confirming F_ = for 6 = 7r/4, as expected, where Fq is defined as the single dot 

spontaneous emission rate. We note that F_ decreases as 9 becomes smaller, suggesting 
that keeping small during the gating operation should be advantageous. However, as 
it is the energetic difference between the unperturbed |00) state and the |— ) states in 
Til, Till, and 7^2 that leads to the CPHASE operation, a smaller mixing angle entails 
prolonging the gating time. The interesting case to consider is that of weak mixing, 
G <C 7r/2, realized by keeping Q/A small. Rewriting Eq. (fT2!) as 

r^ = ^|i-^ ^ I (13) 




{V2Q/A)\ 

and expanding in a Taylor series around Q/A = 0, yields 
Fo fftV „ fft 



For fixed Q we see that F_ ~ A"^ to leading order. 

The CPHASE gate requires that the phase accumulated on |11) exceeds that on 
|01) and |10) by vr, the build-up being due essentially to Vp. During the pulse, |11) 
follows with an energy well approximated by Et^_ = (A — VA^ + 2Q'^)/2, whereas 
|01) and 1 10) follow the dressed eigenstates in their respective subspaces, each of which 
has an energy En = {A — Vp — a/(A — Vp)^ + r2^)/2. Again, in the limit of a large 
detuning, we Taylor expand E(^_ around Q/A = 0, yielding 
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Similarly, expanding around Q/ (A — Vp) = gives 



Therefore, assuming Q is fixed for the moment and that phase accumulates as = 6E t, 
we obtain to leading order 

SE = E,_-2E,^-^J^^, (17) 

for the relevant energy shift during the CPHASE operation. For A ^ the 
denominator is clearly dominated by A^, such that the required gating time t would 
simply take the form t = 7t/{E(^_ — 2E£,) ~ A^/Vf, in the case of a square pulse. Of 
course, the adiabatic gate relies on a Gaussian pulse profile, though in the limit A 3> 
we would expect the above analysis to provide a good approximation. Further, the 
numerical data in the inset of Fig. H] illustrate that the Gaussian pulse duration r does 
also follow a parabolic form with A. Comparing with r„ ~ A~^ from before, we see 
that any decrease in the decay rate r_ by means of a larger detuning A is compensated 
for by the longer duration of the operation. Hence, contrary to the naive expectation, 
in this parameter regime applying a farther-detuned laser pulse does not improve the 
robustness of the adiabatic CPHASE gate to spontaneous emission. 

To validate this result we perform a numerical simulation of the laser-driven 
CPHASE operation [as described by Eqs. (IE]) and (111)] accounting for radiative decay. 
The appropriate ME describing spontaneous emission from a 2LS dressed by a laser 
mode is [16] 

where H = Hs + Hl + is the joint Hamiltonian of the QD (Hs), the laser mode 
{Hl), and the interaction between them {V). Here, x is the joint QD-laser density 
matrix and cr+ and (j_ are, respectively, the raising and lowering operators of the 2LS. 
Any thermally induced emission and absorption processes are neglected as the bosonic 
thermal occupancy N{ujo) is extremely small even at room temperature. Information 
about the state of the laser mode is no longer needed, allowing us to sum over all 
manifolds and obtain a ME for the density matrix p of the system alone. For the 
coherent part of the dynamics, we again retrieve the same form of H2hs as the qubit 
Hamiltonian Hs [16]. The incoherent, spontaneous emission part is described by the 
operator cr_ with associated rate F [IHl HS] . For p the ME can be written in a compact 
way as [50] 

p = -i[Hnh, p] + LsppLlp, (19) 

where iJnh = Hs — {i/2)L\pLsp is an effective non-Hermitian Hamiltonian, and 
is a Lindblad operator for the spontaneous emission processes. Radiative 
recombination affects each dot separately in Tii and TYi', and the Lindblad operator 
describing spontaneous emission in both subspaces is therefore the same as it would 
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Figure 4. Purity of the density matrix for a CPHASE gate witli Tq = 1/100 ps~^, 
Vxx = 5 meV and Vp = +0.85 meV. The dynamic gate uses il = 0.1 meV, whereas 
the adiabatic gates are performed at fio = 1 meV with differing A (given in meV). 
The purity of the adiabatic operation does not depend on the detuning despite the 
substantial differences in the pulse duration, which are shown in the inset. This is due 
to the clearly discernible quadratic dependence of r, the Gaussian pulse duration, on 
A (see text). 

be in the single dot case, L = vT^|01)(0X| and L = ^/T^\10){X0\, respectively. In 
the larger subspace Ti.2 transformation to the diagonal basis still gives a single effective 
Lindblad operator, L = (^/2T qY ^ {ip^\, acting with respect to the Hamiltonian of 
Eq. (gDJ 

We now use Eq. ( fT9il to simulate the adiabatic CPHASE gate at the same time 
comparing it with the dynamic operation proposed in Ref [T6j. We present our results 
as the purity of the full 9LS density matrix, including all decoupled subsystems, taking 
a general input state of the form |0) = (|00) + |01) + |10) + |11)) /2. The purity is 
defined as the trace of the squared density matrix tr[p^]. By definition, it is equal to 
one as long as the system is in a pure state, and decays to l/rf for a totally mixed state, 
where d is the dimension of the Hilbert space [1]. Decoherence effects render the system 
into a mixed state, thereby indicating corruption of the output through a reduction of 
the purity. Figure H] shows the results obtained for a natural excitonic lifetime of 0.1ns. 
The adiabatic curves are plotted for f2o = 1 meV and for various A. We note that they 
nearly coincide, clearly showing that the purity has no apparent dependence on the 
value of A, as expected from our earlier simple analysis. We see also that the adiabatic 
gate retains a purity much closer to unity than the dynamic gate, the speed of which is 
limited by the conditions set out in Ref. [16]. 

I No differentiation between positive and negative values of Vp is necessary throughout this section as 
the radiative decoherence channel is insensitive to the sign. 
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4. Interaction with phonons 

In contrast to an isolated atom, a QD is embedded in the macroscopic solid state 
matrix of its substrate. The equilibrium positions of the lattice ions are determined by 
the overall electrostatic potential of all charge carriers. As a consequence, the electronic 
state of the QD couples to the lattice vibrations (phonons) [51]. Three processes result 
from the interaction of an exciton with phonons in semiconductor QDs: pure dephasing, 
phonon emission, and phonon absorption. 

For many systems, pure dephasing is dominant at very low temperatures j52 | [53 | [M]. 
On excitation, an abrupt change in the charge configuration shifts the equilibrium lattice 
positions of surrounding ions and a phonon packet is emitted as a lattice relaxation 
occurs. This corresponds to a transfer of information about the electronic state of the 
QD into the environment and can be interpreted as a kind of environmental 'which 
way' measurement of the charge qubit [55], resulting in a loss of coherence and forcing 
the system into a mixed state. Pure dephasing is energy conserving and does not 
lead to a relaxation in the QD. It is also an intrinsically non-Markovian process [11] 
and its influence has been studied extensively for Rabi oscillations [56l[57|, absorption 
line shapes [5^ [53] and in the context of the spin-boson model [52] . Several authors 
have shown how pure dephasing may be eliminated if the coherent system excitation 
is slow enough such that all lattice ions adiabatically follow their equilibrium positions 



On the other hand, phonon absorption and emission correspond to real transitions 
between energy levels in the QD. Absorption is only possible at finite temperatures but 
can become a substantial problem in semiconductor structures even for temperatures 
well below lOK. Much like spontaneous photon emission, phonon-mediated relaxation 
processes are possible at any temperature. To study these effects, we derive a Markovian 
ME from first principles. Since we are interested primarily in the limit of low frequency 
dynamics we neglect optical phonons, which are separated by a large energy gap of 
30 meV or more [54]. In this situation, deformation and piezoelectric coupling to 
longitudinal acoustic (LA) phonons are the dominant phonon decoherence mechanisms. 
The interaction between charge carriers and acoustic phonons is generically given by a 
sum over the phonon modes q [51] 



where is the mass density, V the lattice volume and cUq the phonon frequency, with 
wavevector q. The coupling constant is denoted Cq = Dq for the deformation potential 



[13 [m El EH Eg 




(20) 



q 



Here, Mq is the coupling element. 




(21) 
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and Cq = P for piezoelectric coupling, while g is the charge density operator 

^(q) = J2 j dVe-^'i-^^J(r)V;.(r), (22) 



with d\^di the creation and annihilation operators of charge carrier and il)i{r) the 
corresponding wavefunction. For strongly bound excitons with no wavefunction overlap 
between carriers on adjacent QDs, the coupling elements are diagonal and given by the 
difference of those for electron and hole, Mq and Mq, respectively, multiplied by the 
Fourier transform of the electron and hole density operator P['?/'^(r)] and V[^^{y)\. 

= 5^(M^P[<(r)] - MiV[i;1ir)])cU. (23) 

i 

Here, c| (q) denote excitonic creation (annihilation) operators. 

Consider the two-dot system Hamiltonian of Eq. (jl]), once more neglecting the 
biexcitonic level |XX) [§| 

H; = iA- 2\/^)|^_)(^_| + A|^+)(^+| + + h.c). (24) 

The full Hamiltonian including a common bath of phonons is given by if = HI^+Hb+Hi, 
where the bath and interaction terms are given by 

Hb = ^^ CUqa^Oq, (25) 

q 

and 

Hi = Y. 4xCox + 4oCxo) (aq + alq), (26) 
q 

respectively. In the absence of \XX) the phonon interaction term involves only single 
exciton levels, with g]^ = M^p^{ipi) — M^p^{ipi). If the wavefunctions tpQxij^) and V'xo(r) 
are of identical form, albeit centred at different positions ±d, we obtain, by the shift 
property of the Fourier transform, 

^ ^+.q.d (M^P[^^(r)] - M^V[^\v)]) , (27) 
gxo ^ ^-idd (M^P[V'"(r)] - M^V[ij\r)]) . (28) 

We now diagonalize Eq. fl2^ and write Hj in the resulting basis {|C--); |C+)) IV^-)}? with 
respective eigenenergies (A - V2Q^ + A^) /2, (A + V^Q^ + A^) /2 and A - 2Vp. Moving 
to the interaction picture with respect to HI^ + Hb (see Appendix), the system operators 
are ordered by their frequencies, giving 



'G{0,A,T,H} 

X 5^(^?f ± g^')ia^e~'--' + al^e'^-'), (29) 
q 

§ In fact, on resonance, this assumption requires \Vxx — "^Vpl 3> \^\. However, for a larger detuning, 
can be larger as well, and we have checked numerically that excitations to the \XX) level remain 
well below 10~^ for all parameters used in the following. 
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where ± = + for u' G {0, A} and — for u' G {T, S}. Here, 

Po = ^(cos^e|C+)(C+| +sin2e|C-)(C-l + (30) 

Pa= -^sin2e|C-)(C+l, (31) 

Px= cose|7A-)(C+|, (32) 

Ps = -sine|C4(V'-|, (33) 
with frequencies 

A = V2fi2 + A2, (34) 

T = 2Vf- ^(A - (V2fi2 + A2), (35) 

S = ^ ( A + ( vWTA2) - 2Vf. (36) 



In this form, the phonon operators Eqns. ( l32l [33l) are only appropriate for the detuned 
adiabatic gate with A > 2 Otherwise, care must be taken to ensure all frequencies 
are greater than zero. Around resonance and for Vp > 0, S switches sign (see Fig. [2]) 
and we need to use Pj with associated frequency |H| instead. On the other hand, for 
Vp < T switches sign requiring the redefinition of Pt to P\- with associated 

frequency |T|. 

A master equation for the system evolution is now derived in the usual way [l8l [59] 
by integrating the von-Neumann equation for the joint density matrix R of the system 
and bath and tracing over the phonon modes. This results in a reduced system density 
matrix obeying 

t 

p = -J dt'tiph {[Hi{t), [Hj{t'),R{t')]]^ . (37) 



The Born-Markov approximation is now performed, which relies on two assumptions. 
First, we assume there is no back-action from the small system on the much larger bath, 
meaning the joint density matrix can be written as a product at all times R = p® ps- 
Second, the bath has no memory and we therefore replace p(t') by p(t), justified for 
rapid bath relaxation [IHl |59] . If the system dynamics occurs on a timescale much faster 
than relaxation due to interactions with the bath, we can perform a RWA [jj] to arrive at 
an interaction picture ME in Lindblad form [IHl [60] 

p = J+(A) ([A^(A) + 1]Z^[Pa]p + N{K)D[Pi]p) 

+ J^{^'){[N{uj') + l]D[P^,]p + N{uj')D[Pl]p), (38) 

..'e{T=} 

with D[L\p = LpL'^ — l/2{L'^Lp + pL^L) usually referred to as the 'dissipator' of the 
ME. Here, N{uj) = {exp{uj/kBT) — 1)^^ describes the thermal occupation of the phonon 

II The RWA requires A,uj' ^ >/±(A), J±{uj') the diagonal Lindblad-type ME only being strictly valid 
in this limit. Fortunately, the RWA assumption is indeed justified when CIq < 1 and A = as in the 
dynamic case or A > We//i s^s will be used in the adiabatic case. 
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Figure 5. Normalised spectral densities J±{uj): left, for deformation potential, and 
right, for piezoelectric coupling (note the different scales for the y axes of the two 
plots). 



Electron deformation potential 14.6 eV 

Hole deformation potential Dh 4.8 eV 

Piezoelectric coupling constant 1.45 eV / nm 

Effective electron mass nie 0.067 mo 

Effective hole mass 0.34 niQ 

Mass density fi 5.3 g / cm^ 

Velocity of sound Cg 4.8 x 10^ cm / s 

Table 1. Material parameters for GaAs [54 1 152 1 160]. 



modes and J±{w) are the spectral densities given by 

Mu) = 27rJ2 + Cr ^(^ - (39) 
q 

.Uco) = 2^Y1 \9f - Cr - ^.)- (40) 
q 

Note that the Lindblad operator Pq has been dropped since the spectral density vanishes 
in the limit oi u = 0. Consequently, such a ME does not describe non-Markovian 
pure dephasing [521 [51] effects in the dressed basis, on account of the Born-Markov 
approximation, but does provide a suitable description of phonon-assisted transitions. 
In the bare basis, the Lindblad operators Pi^' and Pa consist of emission, and absorption 
terms and Pa features an additional dephasing term. 

For the calculation of J±{uj) we choose the simple case of an isotropic harmonic 
confinement potential for each dot. The excitonic ground-state therefore has a Gaussian 
wavef unction [41] 

/ 1 \^ 

^e/.(r)= — ^ e (41) 
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with d^/h = {h/ y/ff\JhcY^'^ giving the width of the wavefunction envelopes, which depend 
on the different effective masses me/h for electrons and holes. The confinement strength 
c is chosen to be c = 8.3 x 10~^J/m^ such that electrons and holes are subject to a 
harmonic potential V = cr"^ with strength 162 meV at r = 2.5 nm from the centre of 
the dot. Assuming a linear and isotropic phonon dispersion cjq = |q| we then obtain 
for the deformation potential 



J±[^) = — 1 ± sine — 

A-KjJiclR^ \LUpJ \ iUp 

where lOy = Cg/R. Both spectral densities are obviously of superohmic form [M], with 
high-frequency cutoff terms uo^^h = \/^Cs/de^h and uOeh = '^c^j ^/dl + d1 related to the 
finite QD size. These terms filter out phonons with wavelengths too short to interact 
with the dots. For piezoelectric coupling we get ohmic spectral densities given by 



J±{^) = — 1 =t sine 

ATTjJCgRuJp \ UJp 



-(w/oJe)^ _ 2g-('^/'^eh)^ _|_ g-(l^M)2 j _ J^^g^ 



All four spectral functions are plotted in Fig. [5] for the parameters given in Table [H and 
with a dot centre-to-centre distance oi R = 7 nm, fulfilling the assumption of negligible 
wavefunction overlap. Crucially, the dots are still close enough to sustain a significant 
Forster interaction and dipolar shift [H]. From Fig. [5] we identify deformation potential 
coupling as the dominant phonon-decoherence mechanism and neglect piezoelectric 
coupling in the following. This agrees with the literature [521 [51] and is primarily due 
to the fact that electrons and holes couple individually and with a different strength 
to phonons through the deformation potential, whereas only the difference in their 
wavefunctions contributes to the piezoelectric coupling. The difference in deformation 
potentials and Dh combined with the slightly different electron and hole cutoff 
frequencies gives rise to the two peaks in the left hand side of Fig. [5] The drop 
down to zero between the peaks is a 'sweet spot' at which point the phonon couplings 
completely cancel out. Typically, self-assembled QDs are disc-shaped with a lateral 
width of 10 — 20 nm and a smaller height of around 5 nm in the growth direction. 
Accounting explicitly for this would introduce an angular dependency into the spectral 
density with the largest value of the cut-off frequency being in the growth direction, and 
having a value similar to that obtained using the isotropic potential model^. Hence, 
the error introduced by the simple model above is likely to be an overestimation of 
the magnitude of the spectral density towards large frequencies, with no considerable 
change in the position of the cutoff. Furthermore, in the limit of very small frequency, 

% Accounting for the anisotropy would also slightly affect the position of the minima of the deformation 
potential spectral density in Fig. [5l which would cease to be identically zero after averaging over the 
angular dependence. 
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Figure 6. The purity of the fuU density matrix during a CPHASE operation at 
T = 5 K. The system is coupled to a phonon bath as described by the ME p8|) . For 
the dynamic operation, we show results for Vp = ±0.85 meV, as indicated by the labels 
"dyn, +Vf" and "dyn, — Vp". The phase build-up of the dynamic gate (shown in the 
inset) should be step-like but is considerably smeared out in the case of the positive 
Vp, indicative of the fact that this approach does not work (sec text). The congruence 
of the terminal value of the adiabatic curve for A = 2 meV and the dynamic curve 
with negative Vp is coincidental. 

the spectral density does not depend on the form of the wavefunctions for electron and 
hole at all because the Fourier transforms then give V[ip^^'^{r)] ^ 1. 

Solving Eq. (l38|) allows us to characterise the performance of both adiabatic and 
dynamic operation of the CPHASE gate in the presence of a phonon bath. We find 
that the dynamic gate, with positive coupling Vp = 0.85 meV, fares poorly even at zero 
temperature due to phonon emission processes. In this case, is the lowest energy 
eigenstate at resonance such that relaxation from both \(_) and into the dark state 

is possible. By the definitions of T and S [Eqs. ( l35i) and ( l36l) ]. it is clear that a 
reduction of both these frequencies simultaneously to values much smaller than 2Vp is 
impossible. The result is that phonon emission processes, being roughly proportional 
to J-i-(2Vf), always operate at a fast rate, eventually transfering all population initially 
in |11) into \^-)- However, if the Forster coupling is negative, with |Vp| > D,/\/2, the 
upwards shift of is sufficiently large to lift it above both |(^_) and at resonance 
and the dynamic scheme then works relatively well at low temperatures. In contrast, 
in the adiabatic approach phonon emission is always prevented for A > 2 |Vf| because 
|(^_) is then the ground state of the system, yielding a perfect gating operation at zero 
temperature. In fact, if Vp is negative this is even true irrespective of its particular value. 
At finite temperature we expect the adiabatic performance to improve with increasing 
detuning due to the rapid reduction in spectral density for large frequencies. 

In Fig. [HI we plot the purity of the full 9LS density matrix at T = 5 K for both 
the dynamic and adiabatic approaches. The purity of the dynamic gate with positive 
Forster coupling drops quickly to a value of around 0.62, corresponding to the expected 
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state where all population has been transferred into |^_). Of course, decoherence also 
occurs in Tii and Tii', and this is included in the calculation as reflected by the slow 
and steady decline of the curve as time progresses. As expected, the performance of 
the adiabatic gate improves dramatically with increasing A by pushing both A and S 
beyond the spectral cutoff 0, meaning J(A) and decay exponentially. Physically, 
this implies that phonons with an energy matching the transitions between the relevant 
system eigenstates no longer interact with the dot because of their short wavelengths. 
This advantageous regime sets in at an energy below 10 meV, allowing us safely to 
assume a linear phonon dispersion and to ignore optical phonon modes. 

5. Landau Zener transitions 

Non-adiabatic LZ transitions between system eigenstates are a further error source in 
the adiabatic scheme and can only be avoided by varying the control parameters more 
slowly. There is, of course, a limit as to how slowly things can be done, since typically 
many operations need to be performed within the lifetime of the electron spin (roughly 
a millisecond [7]). 

For the simplest case of a 2LS with an avoided crossing, the final transition 
amplitude was derived by Landau, Zener and Stiickelberg [191 120], [21] . Unfortunately, 
this theory cannot generally be applied to driven qubits so that obtaining the transition 
amplitude often involves numerically solving the Schrodinger equation. Wubs et al. [22] 
performed a detailed study of LZ transitions in optically controlled qubits, including 
subtleties such as phase effects of the driving laser. They find that if the final transition 
amplitude vanishes after the pulse has finished, i.e. all population returns to the ground 
state, there may still be a significant population transfer at intermediate times. In 
particular, for a symmetric laser pulse and a transition amplitude not much larger 
than one in a thousand at all times during the operation, almost the entire population 
undergoing a LZ transition is transferred back into the |— ) state once the pulse has 
finished. This is illustrated in some of the curves of Fig. [T] This is undesirable as it 
then becomes difficult to accurately predict the final phase achieved during the adiabatic 
operation. Furthermore, any population transferred to the excited state is much more 
susceptible to decoherence processes. 

In general, any time dependent perturbation of a 2LS will lead to a finite transition 
probability between the dressed states. In order to keep this probability sufficiently 
small, we must explore the conditions under which adiabatic following can reasonably 
be expected. An adiabaticity condition for the simple model of a linear detuning sweep 
A = At and constant coupling strength Q/A<^1 has been suggested in Refs. [HI [18]. 
Here, we take a more general approach and assume that both A and Q are time- 
dependent unless explicitly stated. 

To derive an adiabaticity condition for the CPHASE gate we again consider the 

^ T does not increase with larger A, but this turns out to be unimportant as Px couples states which, 
ideally, remain unpopulated. 
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Figure 7. The population remaining in the |— ) (w |C-)) state of the 7^2 subspace 
is shown for different pulse durations r of a laser pulse Q{t) = ^Iq exp[^ (t / t)] and 
i7o = ^ = 1 ps~^. In this case the adiabaticity condition (|49)) requires r ^ ps. 
The final population leakage for the r = 1.8 ps curve is roughly 8.5 percent and off the 
scale of this plot. 
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Figures. Comparison of adiabaticity conditions. Left: the left hand side of inequality 
(|48|) is plotted as a function of detuning and coupling strength. The pulse duration 
T is chosen to generate a phase shift of tt. The blue colour denotes a region where 
LZ transitions are adiabatically suppressed, i.e. where the LHS of ((48|) is much 
smaller than one. Towards the yellow region, it approaches values of 1/100 and the 
approximation starts to break down. The red colour denotes a value of 1/10 or more 
and adiabatic following can no longer be expected. Right: the same for the LHS of 
Eq. ([311) divided by t. We observe the same qualitative behaviour as in the left part 
of the plot, however, this simplified form of the adiabatic condition is more stringent 
than the previous one, confirming that it guarantees adiabatic following, but is not 
necessarily required for it. 
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2LS approximation to Eq.(jl]), valid as before for large positive detuning. In fact, in this 
regime the level spacing between |^_) and the upper dressed state \(x) is significantly 
larger than that between |(^_) and \C+)- Furthermore, the coupling between |(^_) and 
\Cx) originates from the coupling of \ip+) to \XX), but, at large detuning, \(_) only 
contains a small admixture of Therefore, the probability of LZ transitions from 

IC-) to \(x) is doubly small and it is sufficient to consider only those between |^_) and 
\C+)- The transformation to the basis of approximate instantaneous eigenstates (see Eq. 

o, 

^ /cose -sine \ 
\ smB cosB I 

is time-dependent because and A vary with time (recall that 20 = arctan(v^f^/A)), 
while the Hamiltonian transforms as 

H = HU + i (^^Un U. (45) 



Let = 1/2(A ± vAM-2n^) denote the approximate instantaneous eigenenergies of 
the states and \C+)- In this basis 

H = A-|C-)(C-| + A+|C+)(C+I + e (2|C-)(C+I + h.c.) , (46) 

where the off-diagonal terms in H couple the eigenstates and lead to LZ transitions. To 
achieve adiabatic following, the magnitude of these terms must be much smaller than 
the energetic difference between the eigenstates, giving 

|0|<|A+-A-|. (47) 

This is equivalent to the more familiar condition <^ |A+ — A^| [2^. Inserting 

A"*" we arrive at the following adiabaticity condition 

« 1, (48) 



y2(A2 + 2fi2)3/2 

valid for arbitrary field amplitude Q and detuning A. If the temporal evolutions of 
Q and A are known, this condition can be brought into the form -F(Ao, fio) <^ r, 
where F(Ao,fio) is a time-independent relation of known parameters and r gives the 
characteristic time of the applied pulse. For instance, in the particular case considered 
here of a Gaussian field amplitude fl = r2oexp[— (t/r)^] and constant detuning A = Aq, 
an upper bound for adiabaticity is found to be 

v^fio/A^ < r. (49) 

In Fig. [H] we evaluate the adiabaticity conditions of both Eqs. f llS]) and (jlH]), plotting as 
a function of A and Qo- The same conditions should also hold in subspaces Tii and Tiv, 
with the change ^/2^l fi. We see clearly that LZ transitions do not occur in the limit 
A ^ Qq, exactly the same limit as that for which phonon transitions are suppressed, 
allowing us to circumvent both error sources simultaneously. 
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Figure 9. Fidelity of tlie CPHASE gate. Left: fidelity of the dynamic operation for 
a negative Vp = —0.85 meV. The dynamic reference curve gives the fidelity without 
additional decoherence (see text). Right: fidelity of the adiabatic operation for fixed 
r^o = 1 meV. The inset shows the same curves on a different scale to resolve more 
subtle features. 

6. CPHASE gate fidelity 

To bring together tlie results of the previous sections we now calculate the overall fidelity 
of the CPHASE operation, obtained from a numerical solution to a ME that includes 
both spontaneous emission as well as phonon-induced processes (and, of course, allows 
for LZ transitions). We take an input state given by = (|00) + |01) + |10) + |ll))/2, 
which, after the CPHASE operation, should ideally produce an output state = 
(|00) + |01) + |10) — |ll))/2. However, in practice, for both the adiabatic and dynamic 
approaches some phase is picked up in the Hi and Hv subspaces. This needs to be 
"unwound" using two single qubit operations that we do not consider here. We therefore 
define the fidelity as JF = (0'^|p|0^), where \(f)'j) accommodates the additional phases 
on 1 10) and |01), and p is the full density matrix of the system after the gate has 
finished, including all detrimental environmental effects. As before, we use the material 
parameters of Tableland assume a single dot radiative decay rate of Tq = 0.01 ps~^. 

The results are plotted in Fig. [91 The left hand side shows the fidelity of the 
dynamic gate with Vp = —0.85 meV as a function of the coupling strength and for 
different temperatures. The degradation of the fidelity due to detuned Rabi oscillation 
in both Til and Tii' is included and, as expected, the effect is oscillatory and becomes 
more pronounced as Q increases [16] (to further illustrate this point, we have omitted all 
other sources of decoherence for the curve labelled "dyn. ref."). The fidelity is limited 
by the finite excitonic lifetime towards small values of Q, while to the right of the plot 
phonon-induced processes become more important. The best fidelity achieved here is 
roughly 0.95 at absolute zero, and decreases even further at finite temperature. 

The contrasting behaviour of the adiabatic scheme is shown in the rhs of Fig. [91 
Here, the gate fidelity is plotted as a function of the detuning A, for a fixed value of 
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fio = 1 meV. In agreement with the conclusions of the previous sections, the fidehty can 
be substantially increased simply by applying a suitably large detuning to the driving 
laser. In this limit, the fidelity is restricted by spontaneous emission, as the efi^ects 
of LZ and phonon-induced transitions are confined to relatively small values of A. 
Nonetheless, for A > 5 meV the fidelity is greater than 0.985 for all temperatures shown 
and would improve even further for a smaller spontaneous emission rate (the rate used 
here translates into a rather pessimistic excitonic lifetime of 0.1 ns). For interest, the 
inset shows an intermediate peak only visible at higher temperatures, which is related 
to the dip in the spectral density of the deformation potential (see Fig. [5]). 

7. Conclusion 

To summarise, we have performed a realistic decoherence study of both the adiabatic 
and dynamic approaches to exciton-mediated spin manipulation in coupled QDs. 
We have shown that while dynamic gates suffer from rapid decoherence at finite 
temperatures (to an extent that may prohibit fault-tolerant QC), performing off- 
resonant adiabatic manipulations allows us to greatly suppress decoherence during 
the entangling operations. While a trade-off situation to minimise overall decoherence 
arises for the coupled QD dynamic scheme, this is no longer the case in the adiabatic 
scheme. Since no upper bound is imposed on the gate duration by radiative decay, 
a detuned and slow adiabatic operation is suitable to alleviate the adverse effects of 
both phonon-induced decoherence and LZ transitions. In this case, it is the finite 
spin coherence time that sets a bound on the possible gate duration. However, the 
adiabatic gates we consider operate on timescales of around 100 ps, much shorter than 
typical spin coherence times [7j and only about an order of magnitude greater than 
the "fast" dynamic gates. We therefore conclude that adiabatic optical manipulation 
is a remarkably robust method for entangling spin qubits embodied in semiconductor 
nanostructures. 
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Appendix A. Interaction picture transformation 



The system interaction picture operators of Eq. ( !29l) are strictly valid when the system's 
eigenfrequencies are constant in time, as is the case for the dynamic control scheme. This 
Appendix gives a justification for using the same interaction picture operators within 
the adiabatic approach, which has a time- varying system Hamiltonian due to Q(t). We 
start with a qualitative argument followed by a more formal justification. 

The transformation of Eq. (l26l) to the interaction picture as in Eq. (l29ll would 
be strictly correct if the system frequencies A, T, and H were constant. We require 
each of them to change sufficiently slowly in time to achieve adiabatic following. As 
long as At <^ A, etc. A, T, and S can be assumed to be essentially constant at any 
given moment of time, such that a system operator P^i is well described by P^ie^^'* 
as its interaction picture representation (up to a quasiconstant additional phase only 
varying on timescales much larger than 1/uj')). Since the Markovian approximation in 
the derivation of the ME only depends on the current time, it seems only reasonable 
and consistent to use an instantaneous interaction picture operator. 

More formally, we want to transform 

H{t) = Hs{t) + HB + Hj (A.l) 

to the interaction picture with respect to Hs{t) and Hb- Hs{t), Hb and Hj are given 
by Eqs. ( l24l) . ( 125|) and ( 126|) . respectively. The precise interaction picture transformation 
for a constant Hb and a time-dependent Hs{t) is given by the unitary transformation 

U{t) = Texp I -z / drHsir) ] e"^^^*, (A.2) 







where T is the time ordering operator and we have chosen the Schrodinger and the 
interaction picture to coincide at the arbitrarily chosen initial time t = 0. As Hs{t) 
does not commute with itself at different times, we go into the diagonal basis [in a 
similar way as for Eq. P6|) ] and perform the adiabatic approximation, so that the 
system Hamiltonian is simply 

m) = A-(t)|C-)(C-l + A+(t)|C+)(C+l + Ant)l^-)(V^-l, (A.3) 

where the X^'^{t) are functions of time. Hamiltonian (1A.3P now commutes with itself at 
different times and we can drop the time ordering operator from Eq. (1A.2I) . The system 
part Us{t) of the transformation is then 

Us{t) = exp l^-i j drH'sir)^ . (A.4) 

Instead of using Eq. (1A.4P to transform into the interaction picture, we apply the 
following transformation 

Us{t)=exp{-tH'sm (A.5) 
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to obtain transformed system operators P^^/e*"^'* with nearly constant u', which is the 
form we require for the derivation of the Markovian ME. However, since the total 
Hamiltonian transforms as 



rather than the desired H = Hi as the transformed Hamiltonian. The additional con- 



way then fails to capture the effect of the phonon bath on a small part of the system's 
dynamics. Therefore, the error is likely proportional to the overall effect caused by the 
phonon bath rather than being a constant of a given magnitude. Most importantly, the 
approximation made by neglecting the additional contribution to H improves rapidly as 
A increases relative to Vt. In the limit A ^ fi, for which we predict a very high fidelity, 
the approximation is then very good indeed. 

[1] Michael A. Nielsen and Isaac L. Cliuang. Quantum Computation and Quantum Information. 

Cambridge University Press, October 2000. 
[2] Peter W. Shor. Polynomial-time algorithms for prime factorization and discrete logarithms on a 

quantum computer. SIAM Journal on Computing, 26(5):1484~1509, 1997. 
[3] Lov K. Grover. A fast quantum mechanical algorithm for database search. Proceedings, 28th 

Annual ACM Symposium on the Theory of Computing (STOC), 1996. 
[4] Daniel S. Abrams and Seth Lloyd. Simulation of many-body fermi systems on a universal quantum 

computer. Phys. Rev. Lett., 79(13):2586-2589, Sep 1997. 
[5] Daniel Loss and David P. DiVinccnzo. Quantum computation with quantum dots. Phys. Rev. A, 

57(1):120-126, Jan 1998. 

[6] P. Borri, W. Langbein, S. Schneider, U. Woggon, R. L. Sellin, D. Ouyang, and D. Bimberg. 

Ultralong dephasing time in InGaAs quantum dots. Phys. Rev. Lett., 87(15):157401. Sep 2001. 
[7] Miro Kroutvar, Yann Ducommun, Dominik Heiss, Max Bichler. Dieter Schuh. Gerhard Abstreitcr, 

and Jonathan J. Finley. Optically programmable electron spin memory using semiconductor 

quantum dots. Nature, 432(7013):81-84, 2004. 
[8] T. Calarco, A. Datta, P. Fcdichev, E. Pazy, and P. ZoUer. Spin-based all-optical quantum 

computation with quantum dots: Understanding and suppressing decohcrencc. Physical Review 

A (Atomic, Molecular, and Optical Physics), 68(1):012310, 2003. 
[9] A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. DiVincenzo, D. Loss, M. Shermin, A. Small. 

Quantum Information Processing Using Quantum Dot Spins and Cavity QED. Phys. Rev. Lett. 

83(20) :4204-4207, 1999. 

[10] K. Bergmann, H. Theuer, and B. W. Shore. Coherent population transfer among quantum states 

of atoms and molecules. Rev. Mod. Phys., 70(3):1003-1025, Jul 1998. 
[11] K. Roszak, A. Grodecka. P. Machnikowski, and T. Kuhn. Phonon- induced decoherence for a 

quantum-dot spin qubit operated by raman passage. Physical Review B ( Condensed Matter 

and Materials Physics), 71(19):195333, 2005. 
[12] Pochung Chen, C. Piermarocchi, L. J. Sham, D. Gammon, and D. G. Steel. Theory of quantum 

optical control of a single spin in a quantum dot. Physical Review B (Condensed Matter and 

Materials Physics), 69(7):075320, 2004. 




(A.6) 



this gives 




(A.7) 



tribution can be neglected if H'^t <C H'g (this condition is essentially equivalent to the 
condition on the eigenfrequencies, u't <^ uj' , as stated above). The ME derived in this 



Robust adiabatic approach to optical spin entangling in coupled quantum dots 



25 



[13] A. Grodecka, C. Weber, P. Machnikowski, and A. Knorr. Interplay and optimization of dccohcrence 
mechanisms in the optical control of spin quantum bits implemented on a semiconductor 
quantum dot. Physical Review B (Condensed Matter and Materials Ptiysics), 76(20):205305, 
2007. 

[14] X. Caillct and C. Simon. Precision of single-qubit gates based on raman transitions. The European 

Physical Journal D - Atomic, Molecular, Optical and Plasma Physics, 42(2):341-348, 2007. 
[15] E. M. Ganger, S. C. Benjamin, A. Nazir, and B. W. Lovett. Phys. Rev. B, 77:115322, 2008. 
[16] A. Nazir, B. W. Lovett, S. D. Barrett, T. P. Spiller, and G. A. D. Briggs. Selective spin coupling 

through a single exciton. Physical Review Letters, 93(15):150502, 2004. 
[17] Eliana Biolatti, Irene D'Amico, Paolo Zanardi, and Fausto Rossi. Electro-optical properties of 

semiconductor quantum dots: Application to quantum information processing. Phys. Rev. B, 

65(7):075306, Jan 2002. 

[18] Brendon W. Lovett, Ahsan Nazir, Ehoud Pazy, Sean D. Barrett, Timothy P. Spiller, and 
G. Andrew D. Briggs. Quantum computing with spin qubits interacting through delocalizcd 
excitons: Overcoming hole mixing. Physical Review B ( Condensed Matter and Materials 
Physics), 72(11):115324, 2005. 

[19] L.D. Landau. Zur Theorie der Energieiibertragung bei Stossen. Phys Z. Sowjetunion, 2(4:6), 1932. 

[20] Clarence Zener. Non-adiabatic crossing of energy levels. Proceedings of the Royal Society of 
London. Series A, Containing Papers of a Mathematical and Physical Character, 137(833):696- 
702, sep 1932. 

[21] E. C. G. Stiickelberg. Theory of inelastic colhsions between atoms. Helv. Phys. Acta, 5:369-422, 
1932. 

[22] Martijn Wubs, Keiji Saito, Sigmund Kohler, Yosuke Kayanuma, and Peter Hanggi. Landau-Zener 
transitions in qubits controlled by electromagnetic fields. New Journal of Physics, 7:218, 2005. 

[23] T. H. Stievater, Xiaoqin Li, D. G. Steel, D. Gammon, D. S. Katzer, D. Park, C. Piermarocchi, and 
L. J. Sham. Rabi oscillations of excitons in single quantum dots. Phys. Rev. Lett., 87(13):133603, 
Sep 2001. 

[24] A. Zrenner, E. Beham, S. Stufier, F. Findeis, M. Bichler, and G. Abstreiter. Coherent properties 

of a two-level system based on a quantum-dot photodiode. Nature, 418:612-614, August 2002. 
[25] A. J. Ramsay, R. S. Kolodka, F. Bello, P. W. Fry, W. K. Ng, A. Tahraoui, H. Y. Liu, M. Hopkinson, 

D. M. Whittaker, A. M. Fox, and M. S. Skolnick. Coherent response of a quantum dot exciton 

driven by a rectangular spectrum optical pulse. Physical Review B ( Condensed Matter and 

Materials Physics), 75(11):113302, 2007. 
[26] C. R. Kagan, C. B. Murray, M. Nirmal, and M. G. Bawendi. Electronic energy transfer in cdse 

quantum dot solids. Phys. Rev. Lett, 76(9):1517-1520, Feb 1996. 
[27] Andrew J. Berglund, Andrew C. Doherty, and Hideo Mabuchi. Photon statistics and dynamics of 

fluorescence resonance energy transfer. Phys. Rev. Lett., 89(6):068101, Jul 2002. 
[28] S. A. Crooker, J. A. HoUingsworth, S. Tretiak, and V. I. Klimov. Spectrally resolved dynamics 

of energy transfer in quantum-dot assemblies: Towards engineered energy flows in artificial 

materials. Phys. Rev. Lett., 89(18):186802, Oct 2002. 
[29] Xiaoqin Li, Yanwen Wu, Duncan Steel, D. Gammon, T. H. Stievater, D. S. Katzer, D. Park, 

C. Piermarocchi, and L. J. Sham. An All-Optical Quantum Gate in a Semiconductor Quantum 

Dot. Science, 301(5634):809-811, 2003. 
[30] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, K. C. Nowack, T. Mcunier, L. P. 

Kouwenhoven, and L. M. K. Vandersypen. Driven coherent oscillations of a single electron spin 

in a quantum dot. Nature, 442(7104):766-771, 2006. 
[31] J. Berezovsky, M. H. Mikkelsen, N. .G. Stoltz, L. A. Coldrcn, D. D. Awschalom. Picosecond 

Coherent Optical Manipulation of a Single Electron Spin in a Quantum Dot. Science, 

320(5874) :349-352, 2008. 

[32] J. M. Elzerman, R. Hanson, L. H. Willems van Bcvcren, B. Witkamp, L. M. K. Vandersypen, 
and L. P. Kouwenhoven. Single-shot read-out of an individual electron spin in a quantum dot. 



Robust adiabatic approach to optical spin entangling in coupled quantum dots 



26 



Nature, 430(6998) :431-435, 2004. 
[33] A. Grcilich, R. Oulton, E. A. Zhukov, I. A. Yugova, D. R. Yakovlcv, M. Bayer, A. Shabacv, Al. L. 
Efros, I. A. Merkulov, V. Stavarache, D. Reuter, and A. Wieck. Phys. Rev. Lett, 96(227401), 
2006. 

[34] X. Xu, Y. Wu, B. Sun, Q. Huang, J. Cheng, D. G. Steel, A. S. Bracker, D. Gammon, C. Emary, 
and L. J. Sham. Phys. Rev. Lett., 99(097401), 2007. 

[35] M. Bayer, A. Kuthcr, A. Forchel, A. Gorbunov, V. B. Timofeev, F. Schafer, J. P. Reithmaier, 
T. L. Reinecke, and S. N. Walck. Phys. Rev. Lett., 82:1748, 1999. 

[36] Qianghua Xie, Anupam Madhukar, Ping Chen, and Nobuhiko P. Kobayashi. VerticaUy self- 
organized inas quantum box islands on gaas(lOO). Phys. Rev. Lett., 75(13):2542-2545, Sep 
1995. 

[37] E. Pazy, E. Biolatti, T. Calarco, I. D'Amico, P. Zanardi, F. Rossi, and P. ZoUer. Spin-based optical 

quantum computation via pauli blocking in semiconductor quantum dots. EPL (Europhysics 

Letters), 62(2):175-181, 2003. 
[38] Gang Chen, N. H. Bonadco, D. G. Steel, D. Gammon, D. S. Katzer, D. Park, and L. J. Sham. 

Optically Induced Entanglement of Excitons in a Single Quantum Dot. Science, 289(5486):1906- 

1909, 2000. 

[39] T. Yokoi, S. Adachi, H. Sasakura, S. Muto, H. Z. Song, T. Usuki, and S. Hirose. Polarization- 
dependent shift in excitonic Zeeman spHtting of self-assembled Ino.75Alo.25As / Alo.3Gao.7As 
quantum dots. Physical Review B (Condensed Matter and Materials Physics), 71(4):041307, 
2005. 

[40] T Forster. Discuss. Faraday Soc, 27:7, 1959. 

[41] Ahsan Nazir, Brendon W. Lovett, Sean D. Barrett, John H. Reina, and G. Andrew D. Briggs. 

Anticrossings in Forster coupled quantum dots. Physical Review B ( Condensed Matter and 

Materials Physics), 71(4):045334, 2005. 
[42] Brendon W. Lovett, John H. Reina, Ahsan Nazir, and G. Andrew D. Briggs. Optical schemes for 

quantum computation in quantum dot molecules. Phys. Rev. B, 68(20) :205319, Nov 2003. 
[43] Avinash KoUi, Brendon W. Lovett, Simon C. Benjamin, and Thomas M. Stace. All-optical 

measurement-based quantum-information processing in quantum dots. Phys. Rev. Lett., 

97(25):250504, Dec 2006. 

[44] T. Unold, K. Mueller, C. Lienau, T. Elsaesser, and A. D. Wieck. Phys. Rev. Lett., 94:137404, 
2005. 

[45] Robert Alicki, Michal Horodecki, Pawel Horodecki, Ryszard Horodecki, Lucjan Jacak, and Pawel 
Machnikowski. Optimal strategy for a single-qubit gate and the trade-off between opposite types 
of decoherence. Physical Review A (Atomic, Molecular, and Optical Physics), 70(1):010501, 
2004. 

[46] Claude Cohen- Tannoudji, Jacques Dupont-Roc, and Gilbert Grynberg. Atom — Photon 
Interactions: Basic Processes and Applications. Wiley-Interscience, 1 edition, 3 1992. 

[47] B. R. MoUow. Power spectrum of light scattered by two-level systems. Phys. Rev., 188(5):1969- 
1975, Dec 1969. 

[48] H. P. Breuer and F. Petruccione. The Theory of Open Quantum Systems. Oxford University, 
2002. 

[49] A. Kowalewska-Kudlaszyk and R. Tanas. Generalized master equation for a two-level atom in a 
strong field and tailored reservoirs. Journal of Modern Optics, 48:347-370, February 2001. 

[50] Ulrich Hohenester. Optical properties of semiconductor nanostructures: decoherence versus 
quantum control, 2004. 

[51] Gerald D. Mahan. Many Particle Physics (Physics of Solids and Liquids). Springer, 3 edition, 1 
2000. 

[52] Ehoud Pazy. Calculation of pure dephasing for excitons in quantum dots. Semiconductor Science 

and Technology, 17:1172, 2002. 
[53] T. Takagahara. Theory of exciton dephasing in semiconductor quantum dots. Phys. Rev. B, 



Robust adiabatic approach to optical spin entangling in coupled quantum dots 27 

60(4):2638-2652, Jul 1999. 
[54] B. Krummhcucr. V. M. Axt, and T. Kuhn. Theory of pure dcphasing and the resulting absorption 

line shape in semiconductor quantum dots. Phys. Rev. B, 65(19):195313; May 2002. 
[55] K Roszak, P Machnikowski, and L Jacak. Which way interpretation of the dcphasing of charge 

qubits in quantum dots. Journal of Physics: Conference Series, 30:52-55, 2006. 
[56] J. Forstner, C. Weber, J. Danckwerts, and A. Knorr. Phonon-assisted damping of rabi oscillations 

in semiconductor quantum dots. Physical Review Letters, 91(12):127401, 2003. 
[57] Pawel Machnikowski and Lucjan Jacak. Damping of rabi oscillations in quantum dots due to 

lattice dynamics. Semiconductor Science and Technology, 19(4):S299-S300, 2004. 
[58] P Machnikowski and L Jacak. Dynamical dcphasing of optically controlled charge qubits in 

quantum dots. Journal of Physics: Conference Series, 30:25-29, 2006. 
[59] Crispin Gardiner and Peter ZoUer. Quantum Noise. Springer Series in Synergetics. Springer, 

Berhn, 2nd edition, 2000. 

[60] T. M. Stace, A. C. Doherty, and S. D. Barrett. Population inversion of a driven two-level system 
in a structureless bath. Physical Review Letters, 95(10):106801, 2005. 

[61] A. J. Leggett, S. Chakravarty, A. T. Dorsey, Matthew P. A. Fisher, Anupam Garg, and W. Zwerger. 
Dynamics of the dissipative two-state system. Rev. Mod. Phys., 59(l):l-85, Jan 1987. 



